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AEROMECHANICAL STABILITY OF HELICOPTERS WITH A BEARINGLESS 


MAIN ROTOR - PART I. EQUATIONS OF MOTION 

Dewey H. Hodges* 

Ames Research Center 
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U.S. Army Aviation R&D Command 


SUMMARY 


Equations of motion for a coupled rotor-body system are derived for the 
purpose of studying air and ground resonance characteristics of helicopters 
that have bearingless main rotors. For the fuselage, only four rigid body 
degrees of freedom are considered; longitudinal and lateral translations, 
pitch, and roll. The rotor is assumed to consist of three or more rigid 
blades. Each blade is joined to the hub by means of a flexible beam segment 
(flexbeam or strap). Pitch change is accomplished by twisting the flexbeam 
with the pitch-control system, the characteristics of which are variable. 

Thus, the analysis is capable of implicitly treating aeroelastic couplings 
generated by the flexbeam elastic deflections, the pitch-control system, and 
the angular offsets of the blade and flexbeam. The linearized equations are 
written m the nonrotating system retaining only the cyclic rotor modes; thus 
they comprise a system of homogeneous ordinary differential equations with 
constant coefficients. All contributions to the linearized perturbation equa- 
tions from inertia, gravity, quasi-steady aerodynamics, and the flexbeam 
equilibrium deflections are retained exactly. Part II describes a computer 
program based on these equations of motion. 


1. INTRODUCTION 


The general problem of helicopter aeroelastic stability involves coupling 
between the motion of the individual blades through control system dynamics 
and the rotor wake, as well as coupling between the rotor and fuselage of the 
helicopter The complexity of the problem poses a challenge to the analyst, 
both in developing an analytical model and in understanding its physical 
behavior. An important part of analyzing the general rotor-body dynamic sys- 
tem involves the study of the dynamic behavior of rigid-body fuselage motions 
coupled with the rotor motion. The well-known phenomena of air and ground 
resonance occur m this way. The low-frequency rotor modes interact with the 
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rigid-body fuselage motions to produce xnstabxlxtxes both xn the axr and on 
the ground (refs. 1-4). 

Helicopter rotors without hinges to allow rigid-body flap and lead-lag 
motions of the blades are commonly called hingeless rotors. Hingeless rotors 
have distinct advantages over the more common hinged (articulated) rotors; 
these include fewer moving parts, light weight, and more control— power . In 
recent years there have been efforts to develop hingeless rotors without pitch 
bearings to further reduce mechanical complexity and weight. These rotors, 
called bearingless rotors, rely on a torsionally soft portion of the blade 
called the flexbeam (or strap) that is twisted by the pitch-control system to 
provide changes in pitch. There have been several totally different designs 
proposed within the helicopter industry (refs. 5-9). Of these designs, one 
has a snubber (ref. 6); some have cantilevered pitch-arms without snubbers 
(refs. 5, 7, 8); and one has a torsionally stiff torque tube that is very 
flexible in bending and acts much like a speedometer cable (ref. 9). Both 
hingeless and bearingless rotors have had some aeroelastic stability problems . 
In fact, most production hingeless rotors still must rely on, an auxiliary 
lead-lag damper to suppress ground and air resonance instabilities One 
reason for this is believed to be the lack of suitable analytical capability. 

The analytical treatment of hmgeless rotor air and ground resonance has, 
for the most part, been limited to an equivalent-hinge, spring-restrained 
rigid blade model for the rotor blades (refs. 1-4). The need arises when 
applying these analytical models to somehow arrive at the proper orientations 
for the spring-restrained hinges so as to give the proper values for the aero- 
elastic couplings that arise due to blade elastic deflections, angular offsets 
such as precone or sweep, and the pitch-control geometry. Even with this 
apparent drawback, these relatively simple analytical models have been success- 
ful, to some degree, in assessing the air and ground resonance characteristics 
of some configurations. Bearingless rotors, however, have coupling parameters 
that change significantly as a function of the operating condition. For 
example, when the rotor is at high thrust the flexbeam may be highly twisted; 
whereas, at low thrust it may be untwisted. An alternative to a rigid-blade 
formulation is that of an elastic blade. Again, there is the need to do a 
separate analysis, here to obtain the free vibration modeshapes. This can be 
time consuming for a bearingless rotor because the coupled modeshapes will 
change as a function of operating condition. Also, these analyses are much 
more complicated than the rigid-blade approximations Bielawa has developed 
such an analysis for the bearingless rotor, but only for the hub-fixed case 
(ref. 10). Johnson has developed a rotorcraft aeroelastic stability program 
(ref. 11) that possibly could be adapted to treat bearmgless rotors; however, 
this modification has not been made to date An analysis that possesses the 
simplicity of the rigid blade model but integrates the treatment of aeroelastic 
couplings would appear to be very useful. Such an analysis has been developed 
and is the subject of this report. 

In this report, equations of motion are derived that are suitable for use 
m studying both air and ground resonance as well as the hub-fixed dynamic 
behavior of helicopters with a bearingless mam rotor. The blade is modeled 
as rigid but connected to the hub with a structurally flexible appendage to 
simulate the flexbeam portion of a bearingless rotor blade. The analysis is 
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restricted to rigid -body fuselage motions and all forces that would act along 
the flexbeam are neglected. Only rotor cyclic modes and body pitch, roll, and 
horizontal translations are considered. Each blade has six degrees of freedom 
m the rotating reference frame. This means a total of sixteen degrees of 
freedom m the analytical model. The analysis is tailored to treat the con- 
figurations mentioned in references 5-9 m an approximate way 

In section 2, the physical model is described. In the text, the config- 
uration described in reference 9 is analyzed; modifications necessary to treat 
other configurations are given m the appendixes. In section 3, Lagrange's 
form of D'Alembert's principle is introduced as a means of deriving the equa- 
tions of motion. The coordinate systems used m the derivation are defined in 
section 4 and expressions for certain kmematical quantities are derived. The 
generalized active and inertia forces are derived m sections 5 and 6. The 
equations of motion are formed, linearized about equilibrium, and then written 
m the fixed system in section 7. Particular attention is given to the flex- 
beam structural representation in this analysis. In the appendixes are found 
the details of including blade and body structural damping and the modifica- 
tions to the analysis that are required to treat the different pitch-control 
geometries. 

An independent derivation of the equations of motion was performed by 
Dr. Donald L. Kunz of the Aeromechanics Laboratory m order to check the pres- 
ent derivation. His effort is gratefully acknowledged. 


2. PHYSICAL MODEL 


In this section the physical model used to represent the helicopter is 
described. Only those elements believed necessary to model air and ground 
resonance phenomena are retained. The aircraft dynamical system is composed 
of two parts: the fuselage and the rotor. The fuselage is assumed to be a 

rigid body. When m contact with the ground, the fuselage is suspended by a 
spring system to simulate the elastic restraint to fuselage motion m an 
actual helicopter imposed by the landing gear system. When the aircraft is 
airborne in hovering flight, it is unrestrained elastically. The rotor con- 
sists of three or more rigid blades attached to the hub by means of slender 
elastic beam segments. Both the fuselage and rotor are described in more 
detail below. 

A schematic of the fuselage/hub is shown m figure 1. The hub, mast, and 
landing gear are all included in the mass and inertias of the fuselage. The 
total fuselage mass is mf and the moments of inertia for the mass center are 
I x and Iy, respectively, for the X and Y directions. The aircraft reference 
center, shown in figure 1, is a distance z above the body mass center and a 
distance h below the hub center. In the study of air resonance (in hover) 
and ground resonance, vertical translation and yaw rotation of the body 
uncouple from the other body motions and are, thus, not significant. The 
other four body degrees of freedom X, Y, $ x , and §y are included. These 
quantities are time integrals of velocity and angular velocity components m 
the body fixed axis system, described m more detail m section 4. The 
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landing gear provides stiffness restraining the body motion depending on the 
landing gear geometry shown in figure 2. 

The rotor blades are attached to 'the hub with flexible beam segments 
called flexbeams and rotate at constant angular velocity Q. A schematic of 
one rotor blade is shown m figure 3; here, the details of the pitch control 
system are not shown 

The blade pitch angle may be changed by twisting the flexbeam with any of 
three types of pitch control mechanisms or by using built-m pitch angles 0f 
and 0];, . In this report four cases are considered. Case I is the simplest 
case, with no pitch control system at all. This configuration corresponds to 
an experimental rotor set up to operate only at discrete pitch angles or to an 
operational rotor with a disabled control system The equations of motion are 
developed for this case throughout the text. The modifications for Cases II 
to IV are relatively simple and are given in appendixes A-C, respectively. 

Case II, shown in figure 4, has only a flexible cable-like torque tube tor- 
sionally stiff enough to twist the flexbeam. The cable is assumed to be 
flexible enough m bending to put only a pure twisting moment on the tip of 
the flexbeam. This configuration corresponds to that of reference 9 and also, 
approximately, to the pmned-pinned torque tube of reference 5. Case III 
shown in figure 5 corresponds to configurations with a pitch link and canti- 
levered pitch arm described m references 6, 8, and 9 Pitch change is 
accompanied by bending deflections and vice versa, giving sometimes large 
aeroelastic couplings Case IV, shown m figure 6, is identical to Case III 
except that a snubber, intended to reduce aeroelastic couplings, is added 
The snubber is modeled as an additional nonfunctioning pitch-link/pitch-arm 
assembly. The pitch link is connected to the rotating swashplate and the 
snubber link to a point stationary m the rotating reference frame. All the 
pitch mechanism and snubber mechanism flexibility should be lumped into the 
pitch link and snubber link, respectively. The mass and inertia of these 
components should be included with the blade. Parameters for each configura- 
tion are discussed in the appendixes. 


3. DESCRIPTION OF THE DERIVATION 


The equations of motion are derived through an application of Lagrange’s 
form of D’Alembert’s principle (ref. 12). This principle reduces to Lagrange's 
equations when all the degrees of freedom are generalized coordinates. This 
is the case only for the blade degrees of freedom m this analysis. The body 
degrees of freedom are based on components of velocity and angular velocity in 
the body fixed axis system (fig. 1) and thus are not generalized coordinates. 
They are "quasi-coordinates" and require special considerations, The analyst 
may choose to use the special form of Lagrange's equations for quasi- 
coordinates (ref 13), Lagrange's form of D’Alembert's principle (ref. 12) or 
Newton's laws. 

Kane (ref. 12) has summarized the laws of motion m the form of a single 
equation 
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r = 1, 2, . . . , n 


( 1 ) 


F_ + F * = 0 


where F r are the generalized forces associated with gravity, springs, con- 
tact forces, and aerodynamics. The F r * are the generalized inertia forces 
for the n degrees of freedom. The generalized forces F r are defined for 
the above physical model as 



0V Qk 


+ m b 




r = 1, 2, . 


( 2 ) 


The vectors Fg and Mg are the forces and moments acting on the body at a 
certain point P. The vectors and are the forces and moments acting 
on the kth blade at a certain point Qj^. The velocities and V^ k are 

defined at the point P and in an inertial reference frame. The angular 

velocities co B and u k are written for the body and for the kth blade, 
respectively, also with respect to an inertial reference frame. The degrees 
of freedom q r are u^, v^., w^, 0^ for k=l, 2, ...,b, and 

X, Y, <D X , ® y . 

Similarly the generalized inertia force F r * is defined as 






(3) 


where the inertia forces are 


F b * = -m f a B * 
F k * = -ma k * 


(4) 


Tr ^ 

The accelerations are a B of the body mass center and a K of the kth 
blade mass center. The inertia moments are given according^to Euler’s dynami- 
cal equations, written in detail below. The velocities V B and V k are for 
the body and blade mass centers, respectively. 


In section 5, the generalized active forces due to aerodynamics, gravity, 
body springs, the flexbeam structure, and structural damping are derived. In 
section 6, the generalized inertia forces are derived. Before proceeding with 
the actual derivation, however, it is necessary to describe the coordinate 
systems used m the report and to develop expressions for certain velocities 
and angular velocities needed m the derivation. This is done m section 4 
immediately following. 
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4 . COORDINATE SYSTEMS 


As an inertial reference frame, the space-fixed axes ep e 2 , and e 3 
shown in figure 1 are chosen, e 3 being positive down. The velocity of the 
aircraft reference center is defined to be 


XN a + YN b 


where the unit vectors and Ng are parallel to the body-fixed axes 
and Ng of the N^, Ng, Nq axis system also shown m figure 1. The angular 
velocity of the body is defined to be 

“ B = *x»A + *y»B 


Note that X, Y, $ x , and $y are not geometric distances and angles, they are, 
instead, quasi-coordinates (ref. 13). In this report, the quantities X, Y, 
<3> x , and are assumed to be infinitesimal Thus, 


M= 0 


0 


$ _<f> 

_ y x 


where ep e 9 , and e 3 are unit vectors parallel to ep e 2 , and e 3 , 
respectively. 

The rotating axes N-^, $ k, and N 3 ^ associated with the kth rotor 
blade are also shown m figure 1. The angle i|/k = fit + 2m (k - l)/b. The 
associated unit vectors are related by the following transformation 


~ C K 


\ °*k 


01 K 


-I N 


0 <N 


The kth rotor blade is attached to the hub (at hub radius — e) by means 
of an elastic beam segment. This flexbeam is built m at point 0 (canti- 
levered) along the n^, n 2 ^, n 3 ^ axes, shown m figure 7. The associated 
unit vectors are related by the following transformation constructed from a 
sequence of angular rotations Cf , 3f , and Of 






2 ( = - %% c e f % “ s c f s e f s e f { N 2 ( " £ F V 2 


u f p f s f v f 

-Cfl So S r + So S 


0 f s 3 f s s f + s 0 f s ? f -S e f c s f “ s ^f S Pf C0 f c 3 f c e f l N 3 


6 



The flexbeam of the kth rotor blade is capable of all six beam deforma- 
tions at the tip. The vector of translations is 5 = upn-^ + v k n 2 k + w k n 3 k 
There are also three angular rotations, E k , 3 k , and 0 k , so that the flexbeam 
tip is along the n x ^, ny k , n z ^ axes shown m figure 7. The associated unit 
vectors are related by the following transformation constructed from a 
sequence of angular rotations , B k » © k : 


f n x k 


j n y k 

► = 

(»z k 

The 




C B s 


k C k 


S B 


“ S0 k^k C ?k " C0 k^k C0 k C ?k " S ^k^k S0 k 


k 
s, 

k fa k 




; Ce k s e k s c k + S0 k % 


' s e t c c 


k s k 


- s r s 


^ B k Ce k 


< V> k 


K) 


hi 

r 2 " 

-= [ T k 3«l 

k 

[ k 


k 

l 3 J 


I " 3 J 


(10) 


:h blade is built m at point J along the axes N x k , Ny^, N z ^ 
shown m figure 7. The associated unit vectors are related by the following 
transformation constructed from a sequence of angular rotations £q>» $b» 
and 0^ 


N k 


X 


N k 

» — 

y 


N k 


Z 







“ S0 b%^b " %% 

b Sf3 b% + "VS 


C 0 b C S 


S % Se b S0 b 


’ S0 b C? b " V^b 


C0 b S0 b 

Cg b C0 b 



L k) 


k 


n x 


“x 


n k 

y 

■= EBJ 

_ k 

n y 


k 


k 


n 2 


n z 






(ID 


The N x ^, Nyk, N z ^ axes are the blade-fixed axes. Note that all the blades, 
flexbeams, and angular offsets are identical to each other. This is reflected 
m the absence of the subscript (or superscript) k in the [F] and [B] 
matrices. In the undeformed state [T k ] = [I], the identity matrix. The 
following additional matrices are needed m the derivation below 


[G k ] = [B][T k ] ' 

[H k ] = [G k ] [F] 

[T] - [T 3| ?k=?j g k _g j0k=ie f or k=l,2,...,b 


( 12 ) 


[G] - [B][T] 

[H] = [G] [F] 

i 

Extensive use will be made of these matrices m defining velocities and angu- 
lar velocities needed m the derivation. 


The velocities of several points m the system and the angular velocities 
of the body and kth blade are now developed. The velocity of the aircraft 
reference center is, by definition 
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( 5 ) 


V C = xh a + yh b 

The angular velocities of the body and kth blade are, respectively, 

= ( 6 ) 
“ k = *a»A + *y»B + m 3 + ^3 k ~ M c 4 k *2 k - + 0 k tt x k (13) 


The sequence of blade angular rotations, £ k , S k , and Q k , is shown m figure 8 
from which is determined by inspection. Thus, the velocity of the body 

center of mass is 

V B = (X + zl y )S A + (Y - z$ x )M B (14) 

The position vector of the kth blade center of mass from the aircraft refer- 
ence center is 

r k = hN 3 + eN^ + (i + u k )nj k + v k n 2 k + wjji 3 k + Lx B « x k + cx c H y k (15) 
The velocity of k/' may be established by well-known laws of kinematics 

V k * = £ (16) 

i=l 


where 


Vj k - + (w k + c 3 k ) k a) 2 h 

V 2 k - k V 2 h + (u k + Cl k ) k ( 03 h 

V 3 k = k V 3 h + (v k 4- c 2 k ) V 
\. h = -[(X - hl y ) C(Jj - (Y + 

iC 

+ (Y + h^ x )c^ k ]T l2 
kw i h = - (*x c * k - Vv k )p n + 


(v k + c 2 k ) k m 3 h + u k + c 3 k 

»k + c 3 k ^ u l ti + *k + ^2 k 

(u k + + w k + c 3 k 

^)s^]T i3 + [(X “ M y ) s ^ 

= 1 , 2,3 

+ $ y c ^ k ^ F i 2 + ^13 a 


c ± k = b^b Gk i + cx c G 2i 1 = 3,2,3 

■ t !5L ■ !h! s ac k k • 

C i ■ 9C k h + 3B k 6 k + 30 k \ 


x = 1 , 2,3 


i ( 17 ) 
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It xs necessary to express the angular velocity of the kth blade xn the 


N. 


x 


N„ 


N„ 


axxs system 


“ k = “iV + “ 2 \ k + w 3 k » Z k 


( 18 ) 


where 


N 


i k - 2 <= k V + - w + Wu 

J=1 J 


T1 i k c 0 k B x2 s 0 k B x3 c C k G i2 s £ k G il 


(19) 


Also, the velocxty of the flexbeam txp J xs needed. Thxs xs avaxlable 
by xnspection from equations (17) wxth c^ = 0 


V J = £ k V 1 J n 1 k 


x=l 


( 20 ) 


where 


S J - ’V + % + W - V»3 h 


k-rr J _ k v hi* • k h k h 
v 2 V 2 + + u k m 3 W k U 1 


ky J — k-^7 h 


k, , h .. k h 


3 “ + w k + V*1 - «fc 


to. 


( 21 ) 


Fxnally, xn developing the aerodynamxcs , the velocxtxes of all points 
along the blade N x axis are needed. The blade is moving in space and, in 
addition, there is an induced flow field assumed to remain along the e 3 
axxs. The magnitude of the induced inflow velocity is taken from momentum 
theory (as m ref. 14); for simplicity, it is assumed to vary linearly with 
the radius. 


where 


v^ = fi(e+£ + u + x)(f) 0 < x < L 


( 22 ) 


u k = u + u fc (t) 





0 3/4 ~ H 23 


1 


s S n ( 0 3 / 4 ) 


? 


(23) 
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When the induced inflow is superposed on the blade velocity the result is 

3 

k V'' k k 
V = / , U n 

i i 


i=l 


where 


+ u^ + xG^j + (£w k + xG^g ) ^"0)2^ _ (5-v^ + xG 12 ) to 3 
k = + + xG^ 2 + + xG^^tOj* 1 - (fw^ + ^13) ^ 

U 3 k = + \ + xG k 3 + (£v k + xG k 2 ) k u 1 h - Uu^. + xG^)^^ 

(Y + h$ x - v i $ x )s ijJ F il 


U. 


V 

= - 

(x - 

h$ y ■ 


+ 

[*- 

hl y 


+ 

CO 

H 

> 


£u k 

= a 

+ U k 

+ eF 


= \ 

+ eF 

21 

fw. 

k 

11 

+ eF 

31 


i 7' It 


(Y + h$ x - Vi*x )c «Jf i2 


I = 1,2,3 

\ 

II 


( 24 ) 


(25) 


(26) 


5 DERIVATION OF GENERALIZED ACTIVE FORCES 


In this section, the generalized forces due to aerodynamics, gravity, 
body springs, the flexbeam structure, and structural damping are derived 
Some details are omitted for the sake of brevity. 


5.1 Aerodynamic Loads 

The aerodynamic loads are derived from a quasi-steady version of 
Greenberg's equations (ref. 15) for lift and pitching moment. Also included 
is a quasi-steady profile drag contribution The circulatory lift per unit 
length is given by 

L c = 2 ~- [h + Ve + | (1 - 2x a )e] (27) 
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The noncirculatory lift per unit length is 


L 


NC 


= [h + Ve + Ve + (1 - 4x a )eJ 


(28) 


The pitching moment per unit length is 

.3 


M = - 


pac- 

32 


^2Ve(l - 2x a ) 4- Ve(l - 4x a ) + h(l - 4x a ) 
p sVc^x 

+ (c -|- - 2x a + 4x a 2 jeJ 4 2 — ~ £h + Ve + y (1 - 2x a )eJ 


(29) 


The air density is p and the free-stream velocity is V. The airfoil sec- 
tion is pitched at the angle e with respect to the free-stream as shown m 
figure 9. In accordance with small disturbance airfoil theory, we may set 


u = ^U T 2 + U p 2 = V 

IT s - h - Ve 

r i 

The blade airfoil velocity components are Up and Up, expressed m the 
and direction, respectively. Substitution of equations (30) into 

to (29) yields 


(30) 




"T J [-% + f a - 2x a )s] ' 

£fi[-u P + f tt-4x a)e ] 
pac 2 x UUp 3 r 

2 ^ff r [ U ^ (1 " 8x a + 16x a 

-Upd-4x a )+f (l-f x a 4-fx a 


2 ) 



(31) 


(32) 


Next, the total aerodynamic force per unit length on the blade airfoil 
section is considered The noncirculatory lift is taken to act m the direc- 
tion normal to the chord line as shown m figure 10. The circulatory lift is 
taken normal to the resultant velocity U. An aerodynamic profile drag force 
per unit length D, acting parallel to the velocity resultant, is included, 
based on a constant profile drag coefficient cj . 


D = 


pc, cU^ 


(33) 


The force components and directions as shown in figure 10 give the following 
expressions for S and T 
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where, from figure 9 


S = L n sm a - D cos a 

L* 

T = cos a + L^ c + D sm a 


I 


cos a = 


U 


sm a = - 


U 


(34) 


(35) 


Substitution of equations (31), (33), and (35) into equations (34) yields 


pac 

2 


T = 


pac 

2 


Up 2 - f (1 - 2x a ) Up e - UU T J 
-U P U T + f (1 " 2x a )U T e - ~ Up + 



(1 - 4x a )e 





(36) 


Since c^ /a is small with respect to unity and because the magnitude of 
the aerodynamic pitching moment is small, it is permissible to set U = Up m 
S, T, and M. The aerodynamic force per unit length acting along the blade is 



The aerodynamic pitching moment per unit length is 



(37) 


(38) 


Now, for the kth blade 
q — r TT ^ — — fT — 9-v it h k 

S k 2 l_ Up 2 ^ 2x a^ u P 03 1 

T k - £ r [-W - f «p k + f a - 
K k - - p f- (vA 11 + (f) 2 [ IJ i k,i 


- dD i k2 ] 

2x a )H I k » 1 k + - 4x a )i, 1 k ] 

- 8x a + 16x a 2 )wj k - (1 - 4x a )U p k 


(39) 





/ 


where d = c<j Q /a has been neglected with respect to unity, and e has been 
replaced with the component of angular velocity of the blade m the N x k 
direction The velocity components U p k and U T k are given by 
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Ut 


- X) ^3! 


1=1 

3 


I 


u T k = 2 U i 


k P k 

J 2i 


1=1 


where the velocity components are defined m equations (25) . The 

dynamic force and moment resultants at the point J are given by 


w. 


aero 


S f k N y k + T f k N z k 


(*k) - «*\ k - T »V + S n\ k 

v 7 aero 


where 


k _ 


m k 

T f “ 


( 

t 


S k dx 


T k dx 


m 


m 


■( 

■jf 


xS k dx 


xT k dx 


*-jf 


M k dx 


The generalized forces due to aerodynamics now follow immediately from 
equations (2) and (4) 


( F "k) ' ( F 4ero ‘ 

K) ' (^Lro ’ ” 2 

K) = VVaero ' 113 

W = ( Bk L ' “• 

\ /aero . . 

W = ' ' aero ' 

\ /aero 

(\) * We. ' 


Cv- Up S- 11 -I 

, k 


> k = 1, 2, . . . , b 


(40) 


aero- 


(41) 


(42) 


(43) 
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( F *) ■ E (" k ) 

' 'aero ' 


k=l 

b 


aero 


W - E W 

' n ova 


N 


B 


aero 


k=l 

b 


aero 


[%(< 


aero k=l 

+ 


V 11 + V i! ) ni ‘ 


c ,) k F 3l - *32) - *»k(‘ 


c ^ k F 21 - s ^ k F 22 ;[“l 


*)} 


+ [ £Wk ( C ^ k F U ‘ S \ Fl2 ) " £U k( C \ F31 ~ S \ F 32)] n 2 k 
+ ^k^ c ^ k F 21 ~ s iJ; k F 22^ ~ ^k( c t|> k F ll “ s v k F l 2^J n 3 k | 


W 


aero k=l 

+ 


w 
r (w 


* n a . 

aero «■/ 


aero 


E ( V 1 * ■ V i2 H 

1=1 


[*” k ( s * k F 21 + c (, k F 22 


) - ^( S * k F 31 + V 32 


)]”> k 


[ W k( s * k F 31 + ^ k F 3 2 ) " »k( B * k F ll + c * k F l 2 )] n 2 k 
[ lF k( s * k F n + <^12) - + <=* k F 22)]“3 k | 


W. 


Nv 


(43) 

concluded 


aero 


Equations (43) define the contribution of aerodynamics to the generalized 
forces. 


5.2 Gravity 

In this section the gravitational contribution to the generalized forces 
is expressed. For a system of rigid bodies, gravity applies a force acting 
through the mass center of each component. For the body 
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(44) 


e 3 

and similarly, for the kth blade 

< F k) g = m S e 3 (45) 

The subscript g refers to gravity. The symbol g is the rotor thrust per 
unit mass of the aircraft (^/M) if the aircraft is airborne. When the air- 
craft is m ground contact, the analysis must be slightly modified. This is 
done at the end of section 5,3. 


The generalized forces due to gravity follow immediately from equa- 
tion (2): 


( F 0 g ■ 
K) g - 
Wl ■ 


mg e 3 • n x 


mg e 3 ♦ n 2 


mg e 3 • n 3 


K) ■" ge » ‘Si r*i' 


4 3r 

i=I 


(xl = ” se 3 -Y^t< 


sr 3e k 


(F x ) - Mg a, • N a 


<Vg - Mg e 3 • S B 


K) ' mge3 


• — f 36. 

i=l k 


( ? O g = (mgbh ' B f S2)e 3 * % - mg e 3 • £ j [(^k +c 3 k )( 


c ^ k F 21 S i|i k F 22j 


rk + c 2 XV 91 “V^ 


)] n l k+ [( £u k +c l k )( c ^ F 31- s ^ F 3 2 ) / (46) 


^k +c 3 )( C ^ k F 1 1 “ S ^ k F l 2)J n 2 + 


( Av k + C 2 k )(‘ 


,V 11_ V 12 J 


' /\ 'K it /j \ l\ /\ r k ■t'k 

V U k + c l k )( c ^ k F 21 " s i|) k F 22)J n 3 k | 

= - (mghh - m f gz)e 3 • H A + mg e 3 • g J [(^ k + c 3 k ) (s^ 21 + ^22) 

' '% 

~ ( v k +c 2 k )( s <f. k F 31 + c l ( )k F 32)] n l + [( u k +c l k )( s ^ k F 31 + c l J; k F 32) 

- (w k +c 3 k )(s^F 11 + c^F 12 Ji l2 -f [(v k +c 2 k )^ k r n + c^F 12 J 

- (u k + c 3 k ) (s^F 21 + c^F 22 )j n 3 J 


15 



Equations (46) define the generalized forces due to gravity when the air- 
craft is airborne. If the aircraft is hovering at thrust = weight, g becomes 
the acceleration of gravity g. When the aircraft is on the ground, the body 
forces F x and Fy- need slight modification.' This modification is closely 
connected with the body spring forces treated m the next section. 


5.3 Body Springs 


When the aircraft is m ground contact, a system of springs is introduced 
into the mathematical model to account for restraint of fuselage motion due to 
the landing gear flexibility. These springs, shown in figure 2, are symmetric 
about the aircraft reference center both in the longitudinal and lateral 
directions In deriving the generalized forces, the analysis will be carried 


out in only one plane, the 
total system. 


e^ , e 3 plane and the results extrapolated to the 


Consider the two-dimensional system of figure 11. The velocities of 
certain points are needed m order to derive the generalized forces. The 
reference center velocity is 

V C = XN a + ZN C 


(47) 


and the body angular velocity is 


m B = 




(48) 


The velocity of points Aj and A 2 is, respectively, 

1 


= (X + £ z $ y )N A + (z + ~ * 7 )h ( 
V Az - (x + * z *y)N A + (z - ~ $ y )N 


(49) 


(50) 


The force applied by the springs at each of the attachment points is 
% = “ 2K x (e l + *zV e l “ 2K z( e 3 + T $ y) e 3 

F A 2 = ~ 2K x (e i + W e l " 2K z( e 3 ' ~T $ y) e 3 

where e^ and e 3 are the space-fixed deflections of the points A-^ and A 2 
along the ej^ and e 2 axes. The vertical degree of freedom is retained only to 
aid m modifying the gravity terms m F x and F y for the ground contact case. 
When the aircraft is on the ground, there is a steady vertical deflection of 
the ground springs that varies with rotor thrust and gross weight. Thus, for 
first order m X and $y the thrust and weight Mgg must be included 

for this case. Each of these forces may be assumed to act at the reference 
center for the purpose of calculating spring forces. 


16 


ORIGINAL 
OR POOR Q.UALU 



F B* “ " + Mg 0 e 3 


(51) 


From equation (2) 


F x = _4K x (e i + *zV + 4K z e 3 $ y “ Mg o $ y 

F z = -4K Z e 3 - P + Mg 0 (52) 

F * y = -(4K x £ z 2 + K z £ x 2 )$y _ 4 ^ ^ 

For the equilibrium of the coupled rotor-body system (all degrees of free- 
dom = constant), F z = 0. Therefore, 

4K z e 3 = Mg 0 - ?T (53) 


Also, the deflections along e 3 and e 3 are identical to X and Z for infini- 
tesimal $y. Thus, for X and $y degrees of freedom 

(F x) sp r = “ 4K x (X + My> + (Mg 0 " ^>*y “ M S 0 *y ) 


= -4K x (X + ft z *y) - 


(54) 


( F t y ) “ -f 4 V, 2 + Vx'l'y ~ 

spr 


Similarly, for Y and $ x degrees of freedom 


(F y ) 


spr 



= -4Ky(Y - £ z $ x ) + 

= ~(4Ky£ z 2 + K z £ y 2)$ x + 4Ky£ z Y 


(55) 


Equations (54) and (55) express the generalized forces due to body 
springs. The gravitational terms m F x and F y are replaced by rotor thrust 
terms when the aircraft is m ground contact. This is equivalent to replacing 
g m F x and F y with g* where 

g* = -f- (56) 


when the aircraft is on the ground and replacing g m all other terms with 
8 0* 


5.4 Flexbeam Structure 

The structural loads exerted by the flexbeam are represented symbolically 
by the flexbeam stiffness matrix [K] . The matrix [K] is determined numeri- 
cally by perturbing the equilibrium solution. This operation is described m 
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section 7.2 below. With the stiffness matrix formulation, the generalized 
forces are 


(v) 

' ''flex 


i=l 


K q 

J 


k _ 


- VVWW k 1,2,. . .,b 


( 57 ) 


The fact that all the b flexbeams are identical is reflected m the lack of 
dependence of the [K] matrix on k. 

When S2 = 0 and the flexbeam is undeformed, the matrix [K] is given by 


[K] = 


EA 

a 


0 


0 


0 


0 


0 


12EI C 


0 


~6EI C 

£* 


0 


0 


12EI f 

£ 3 


0 


-6EI f 

£ 2 


0 


4EI C 

£ 


0 


0 


4EI f 

£ 


0 


GJ 

£ 


(58) 


where EA, EI C , Elf, and GJ are the flexbeam axial stiffness, chordwise bend- 
ing stiffness, flapwise bending stiffness, and torsional stiffness, respec- 
tively. When the rotor is at some general operating condition, the matrix [K] 
fully couples the blade equations structurally. 


5.5 Structural Damping 

In the coupled rotor-body system only the blade lead-lag motion and some 
body modes have small enough damping from aerodynamics to be significantly 
affected by a small amount of structural damping. The non-zero generalized 
forces due to structural damping are 
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^ ^k^damp 

-c ? C k 

(^ x )damp 

-c x X 

(^y)damp 

-c y Y 

( F$x )damp ” 

x ' 

^damp = 

-c* $- 
y ■ 

C X , c Y , c $2 



( 59 ) 


terms of damping ratios for both blade and body motion. 


This concludes the development of the generalized active forces. In the 
next section the analysis continues with the generalized inertia forces. 


6. DERIVATION OF GENERALIZED INERTIA FORCES 


In this section, the generalized inertia forces due to body motion and 
blade motion are derived. As m the last section some details are omitted for 
the sake of brevity. 


6.1 Motion of the kth Blade 


Generalized inertia forces that are associated with generalized coordi- 
nates, as the blade degrees of freedom are, may be derived directly from the 
kinetic energy. For the rotor system the kinetic energy may be expressed as 


T = 


m 

2 




(60) 


Once this quantity has been expressed m the system degrees of freedom, the 
generalized forces for the blade degrees of freedom are given by 




VVVWk 


k =. 1, 2, . . . , b 


( 61 ) 


This operation is straightforward and the details are not given here. 
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6.2 Motion of the Body 


Unfortunately, the generalized inertia forces associated with body motion 
cannot be derived f-rom equations (61). This is because the degrees of freedom 
describing body motion are not generalized coordinates. They are, instead, 
quasi-coordinates and require special considerations. The portions of body 
generalized inertia forces that are linear m blade variables are already 
known because of symmetry considerations and equations (61) Hence, only 
terms linear in X, Y, $ x , and $y and their time derivatives need to be 
retained m this section. 


The accelerations of B x and k* used in equation (4) to define the 
inertia forces are needed. These may be derived from equations (5) to (17) 
by standard laws of kinematics 


a B * = (X + z$ y )N A + (Y - z$ x )N b + . . . 


.k* - 


where 


/ 3 \“ 


/ 3 \ ~ 

x - (h + Z doXsK 

N A + 

Y + h + Z 

\ i=l / 


L V 1=1 / -I 


N. 


j(^X 2 + F l2 C lJ>k) 

3 

■ (*y + 2M X ) 2 d i k ( F il^ k - F i 2 ^ k ) 




di k = Zu k + c x k I 

d 2 k = £v k + c 2 k 

d 3 k = £ "k + C 3 k 


(62) 


(63) 


Here, the dots refer to blade terms and nonlinear terms m X, Y, $ x , and 3>y 
and their time derivatives. Also, in equation (3) the inertial moments must 
be written. These may be expressed directly from Euler's dynamical relations 
and the angular velocities in equations (6) 


m b* = - Wa - Vy^B + • • • 

*k* = [^2 k “3 k < 1 2 “ I 3 > " “l k IllN x k 

+ [m 3 k 03 1 k (I 3 - Ij) - “2 k ^2^y k 

+ [a) 1 k 0) 2 k (I 1 - I 2 ) - w 3 k I 3 ]N z k 


(64) 
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Now the body terms of the generalized inertia forces for X, Y, $ x , and 
$y may be written from equations (3), (4), (62), and (64) 


■x 


n a • P 


f b* + X r k* 

k=l 


) + . . 


F v* - *B • ( r B* + X 


5 '■■) 


b 

+ X 

k=i 


F * x t - • F B* + % • (v + X «k*) 

fh + E + X d l k ( F H^ k + F l2 <: * l ) N C 

F 4. y * ’ =% • V + ®B ' f«B* + X «k *) 

-(h + X a> u j» A + X - F l2 s * k ) N . 


b 

+ E 

k=l 


* F k 


(65) 


Equations (65) are the body terms of the generalized inertia forces for the 
body degrees of freedom. In the next section, the equations of motion are 
written from equation (1). 


7. FORMATION OF THE EQUATIONS OF MOTION 


The equations of motion are to be written from equation (1) using all the 
components of the generalized forces m sections 5 and 6. The ultimate pur- 
pose of this analysis is to provide a means of assessing the linear stability 
of small motions about equilibrium. In this section the equilibrium solution 
is discussed separately before the final linearized perturbation equations are 
written. Because from the outset the blade and flexbeam properties have been 
assumed identical for all blades, the equilibrium values of u k , v k , w^, ? k , 
3 k , 9k are identical constants for all k. Thus, the blade degrees of free- 
dom have steady and oscillatory components 

u k = u + u k (t) ; ? k = ? + C k (t) 1 


v k = v + 


B k = B + B k (t) 


( 66 ) 


W k = w + W k (t) 


®k = ® + I 
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Pecall that the equilibrium values of X, Y, $ x . 


and $y are zero so that 


X = X(t) 
Y = Y (t) 




x 


l x (t) 


^y = *y( t) 


( 67 ) 


The tilde refers to an infinitesimal perturbation motion. Collection of terms 
from sections 5 and 6 yields a set of relations like the following 

F + F* = F - F ft) = 0 
ru k “ k U u k ^ 


‘v. 


W, 


+ f: 


* _ 


= F,. 


F (t) = 0 


+ F* = F, 7 - F t , (t) = 0 


' w 


w. 


+ F* « F- - t r (t) = 0 


a k 


*k 


F r + F r = F r - F R (t) = 0 

B k 3 k 6 


+ F 


* _ 


= F 6 - F q (t) = 0 


e. 


F x + F X* = - f = 0 


Fy + F y * = - | F y (t) = 0 


+ F* 


X 


‘X 


f \ (t) = 0 


X 


+ FS 


= - \ % (t) = 0 


( 68 ) 


__Here all perturbation quantities have been linearized m u k , v k , w k , £ k , 
f3 k , 6 k , X, Y, $ x , and 4> y and their first two time derivatives Below, the 
rotor equilibrium solution is obtained by performing an iterative static 
structural analysis for the flexbeam. The flexbeam structural stiffness 
matrix is next obtained from a numerical perturbation of the flexbeam equilib- 
rium position. The perturbation equations, although linearized m all the 
perturbation degrees of freedom, contain terms with periodic coefficients 
(8 V c ij; k ) • In order to most efficiently solve the system of linear, ordinary 

differential equations, the system is transformed to fixed (nonrotating) 
coordinates by the so-called multiblade coordinate transformation (ref. 16). 
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7.1 Equilibrium Solution for the Elexbeam 


The equilibrium generalized forces F u , F v , F w , F^, Fg, and Fq are, 
physically, the components from aerodynamics, gravity, blade inertial loads, 
and flexbeam structural loads, of force and moments acting at the tip of the 
flexbeam. Explicit expressions for these generalized forces are known except 
that expressions for the flexbeam structural loads are unknown. This means 
that all loads external to the flexbeam are known explicitly m terms of the 
deflections of the tip. When the details of the generalized forces are 
carried out from sections 5 and 6, equations (43), (46), and (60), the flex- 
beam tip external force Ft and moment M t are given by 

F t = F u n i + Va + *>3 ) 


3 

- E 

1=1 


M n 
i l 


Mi =- 


c S s £ 


F r + s r F R + — F q 
eg C Z B Cf3 0 


S -So S c 

M~ = — Fr - c r F R + — F fl 

2 C B e C S C B 8 


M 3 = F c 


(69) 


The bars over F u , F v , F w , F^, Fg, Fg indicate that these are-external forces 
and moments, excluding the structural part. These quantities are listed below 
m nondimensional form. Forces are nondimensionalized with respect to I& 0 2 /£ 
and moments by IQ Q 2 . 

F u = inft 2 £[£u + c x - F 13 (J 33 - h)] - mg£F I3 + Z(S f ^G 21 + T f ^G 31 ) 

F v = mS 2 £[£v + c 2 - F 23 (J 33 - h)] - mg£F 23 + + T f G 32 ) 

F w = mjj 2 £[lw + c 3 - F 33 (J 33 - h)] - mg^F 33 + J(S f ^G Z2 + T f ^G 32 ) 


F ? = n z 


6 


3H 


13 


13 a c 


+ mQ* 


1=1 


9c! _ Sc 2 

(£u + Cl ) — + (Iv + c 2 ) — 


+ (£w + c 3 ) 


3c 3 3 ^33 


9C 


3? v 33 


(J„-h) 


a j 


33 


m§ 3? +G 13 M m 0 G 2 3 T m 0 + G 33 S m c 


(70) 
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f s . 52 




3 H i 3 

^x 3 — 3 g *■ m ^ 2 


- i=l- 


9Si _ 9 c 2 

(Jtu + Cj) -TT- + (£v + c 2 ) 


93 


9c 3 3J 33 


+ u»+<= 3 > jg 33- (j, 3 -E) 


9 J 


33 


“ mS -JT “ n l Mra n + T1 2 T m n " ^ 3 S M, 


F 0 = Cl 2 


E^ 


9 H. 


13 


13 96 


+ mSl 2 


1=1 


9 c, 


0 ‘ ‘“0 J "0 
9 c 0 


(£u+c 1 ) "37“+ (Jlv+c 2 ) 


96 


\ ( 70 ) 
con- 
cluded 


9 c- 3 J-,g 

+ (J,w+c 3 ) — — (J33 - h) 


9 J 


- mg 


33 


ae + B ll M m 0 B 21 T m 0 + B 31 S m 0 


where 


m = 


SI = 


ml/ 


‘0 


Z = - 
* L 


C i L 


r _ h 
h " L 


J 33 = h + (£u + c 1 )F 13 + (hr + c 2 )F 23 + (£w + c 3 )F 


33 


§=~V 


J - = T 


q _ I 

f 0 6 


V 2 - dV 2 + 3 (U_. V_ - dU T V T + U 2 - duS 


ry V r p 
0 0 


3 c 


- ~r (1 ~ Zxj.w (Vp + 2u 


< l 

>.(% + 2 \) 


,) 


( 71 ) 
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Tf - - ± 


3c 


Ycx = 


M m„ 


6 

Y5* 


P„X + i ( U Po VT » + V P„»T 0 ) + 3U Po“P 

4 a - 2x a )» lo (v Tq + 2U T Jj 

! (VP 0 + V 1 .) 


Vt, V T + 
p 0 0 


+ V U T I + 3U U T 

p 0 0/ p o ± 0 


T - 

■ L m„ a 


» 1 „ (vt 0 + 2 % 0 )tl - 8x a + 16x a 2) 

V P V T +4 ( U P V T + V p U T ) + 2U p U T 
p 0 0 J \ p o ± o p o 1 0 / p o i c 



L 

QH i3 

"^23^31 ^ 1 3^32 

= f2(-J23 H 21 + ^13 H 22 + X ' U< J )H 23-' 

V P 0 - 5 « H 33 - H 23> 

V T„ * »W H 23 + H 33> 

Jl 3 = S(uF 11 + vF 21 + wF 31 ) 

J 23 = A(5F 12 + vF 22 + wF 32 ) 



sssss* 


> (71) 

| con- 
I eluded 


Each of the force and moment components acting at the tip of the flexbeam are 
known functions of u, v, w, 3, and 0 which are unknown. In order to 
solve for u, v, w, £, 3, and 0 we assume a set of values u^, v t , w t , £ t . 
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g and 0 t . We can then calculate the forces and moments F t and Mt based on 
equations (69) and (70). The force Fr and moment Mr at the flexbeam root 
are then determined from statics 


F 


R 


= F* 


% = M t + r t X F t 

r t = (1 + u t )n x + v t n 2 + w t n 3 


* 


% = M t + (v t F w - w^)^ + [w t F u - (1 + u t )F w ]n 2 + [ (1 + u t )F v - v t F u ]n 3 J 


(72) 


At any station s along the deformed flexbeam length, the force F and 
moment M are 

F | 


H = H r - r X F r 


> (73) 


r = [s + u^(s)]n 1 + v(s)n 2 + w(s)n 3 J 


where u^(s) is the geometric component of u(s) due to bending only. The 
bending moments m the principal axis system for the flexbeam cross section 
, u 2 , and u 3 are 


M • u 2 = Ak 
M • u ^ = BA 


) 


(74) 


where A and B are flap and chord flexbeam bending stiffnesses made dimen- 
sionless by IQ 0 2 £. The curvatures are k and A, made dimensionless by 1/2.. 
The torsion moment is, according to reference 17, 


M • n 


1 



(75) 


where, from purely geometric considerations, it can be shown that curvatures 
and slopes are related by the following" 

£ = ~Kc e + As 0 ) 


dz 

ds 


^0 + Xc 9 
c £ 


( 76 ) 


d6 

ds 



/ 
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and where 


E = ^^ a2 _ + B ’ ) - (for a rectangular cross section) 
5D 

T = F • o 1 


(77) 


The torsion stiffness C is made dimensionless by and the axial 

stiffness D is made dimensionless by The longitudinal strain is 

given by 



(78) 


which is equal to, according to reference 17, 



(79) 


We define u(s) as being 


u (s ) = u b (s) + u s (s) 


(80) 


where u s is due to axial stretching and u^ is a purely geometric deflec- 
tion due to bending only. We then set the strain due to bending equal to 
zero and obtain 


e 



(81) 


In combination with equation (79), we have 



(82) 


Also, from purely geometric 


considerations (see the appendix of ref. 17) 


du b 

ds “ c e c s “ 1 


dv 

ds 


Vc 


dw 

ds ” s 8 


> 


(83) 


Now, with initial conditions u s = ub = v = w= t; = g = 0 for s = 0 
equations (76), (82), and (83) may be integrated numerically, using equa- 
tions (73) to (75) and (77), from s = 0 to s = £ + u s (s). This will yield 


27 


ORIGINAL PAG® W 



a set of values for u, v, w, £, B, and 9 at the tip of the flexbeam which 
will be the same as the set of u t , v t , w t , £ t , Bt» an d 0 t assumed if and 
only if the trail set is the correct set. Since the correct u t , v t , Wt> £.t> 
B^j and 9 t are unknown, the following function may be minimized 

min 

u t .,v t .,w t , = <u ~ u t ) 2 + (v- v t ) 2 + (w- w t ) 2 + (? - ? t ) 2 + (B - B t ) 2 + (9 - 0 t ) 2 

5 £ ,B t > 0 t (84) 

using the modified Levenberg— Marquardt algorithm (ref. 18). The minimum 
SF - 0 is the exact solution for the rotor equilibrium deflections. 


7.2 Flexbeam Stiffness Coefficients 

The_external equilibrium generalized forces at the flexbeam tip F u , F v , 
F w , Fg, F 0 = Q-l, i = 1, 2, . . ., 6 may be used to obtain the flexibility 
coefficients for the flexbeam structure. The Q 1 may be perturbed m succes- 
sion, starting with Qj , by a small number e This will yield six different 
sets of the deflections u, v, w, £, B, 9 = x.,^ i = 1, 2, ...» 6 only 
slightly different from the equilibrium values. Flexibility influence 
coefficients are simply 


x (Q + e) - x (Q ) 

F = -^-2 ( 85 ) 

13 £ 

The smaller value of £ that is chosen, the closer the coefficients will be 
to a pure linear perturbation. If the deflections are calculated to N 
significant digits, £ should be chosen at 10~ N ' 2 to minimize both nonlmeari- 
ties and truncation errors, simultaneously. 

The stiffness matrix for the flexbeam is simply 

[K] = [F" 1 ] (86) 

When Q = 0, the structural stiffness matrix is given by equation (58) . The 
numerical scheme described above gives this stiffness matrix to five signifi- 
cant figures when = 0 and exhibits a fully coupled system structurally, 
for general flight conditions. 


7.3 Transformation to Fixed System Coordinates 

Although the perturbation equations are linearized m u^, u^, w^, 4 ^, 

Bk> 0 k> Y, $ x > ®y an ^ their first two time derivatives, the equations have 
periodic coefficients in the form of sm and cos These terms may be 

eliminated with the multiblade coordinate transformation (ref. 16) when b > 3. 
It is necessary to retain only the rotor cyclic modes since the collective 
modes couple only with vertical translation and yaw rotation, which are also 
decoupled from the retained degrees of freedom m hover. The differential 
collective modes and the warping inodes (for b > 5) are uncoupled from all 
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other degrees of freedom m hover. This reduces the 6b rotor blade degrees 
of freedom to only 12 (two cyclic modes per blade degree of freedom); thus a 
total of 16 degrees of freedom remain for the coupled rotor-body analysis. 

The 12 rotor degrees of freedom, m the order that the matrices below are 
written, are 


= f E \ cos *k 

k=l 
2 b 

“s ■ b £ ”k Sln 
k=l 


2 A ~ . 

v c = b E v k COS *k 

k=l 
? b 

*s = b £ \ sin \ 
k=l 

2 b 

= b £ \ cos % 
k=l 


- f E \ sin * k 

k=l 

5c = | E \ 

k=l 
2 b 

5s - b E C k si” * k 

k=l 
2 b 

= b £ e k cos ^k 
k=l 

2 b 

h - b ,E \ Sln *k 

k=l 


§ c = f t \ cos 


k=l 
b 

6- = £ £ \ sin * k 

k=l K 


2 : 

's b 


(87) 


, page IS 
QUALITY 
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The use of the varxables xn equatxons (87) wxll elxmxnate all terms wxth 
perxodxc coeffxcxents xn the body equatxons. The transformation is carrxed 
out on each kth blade equatxon xn the following manner* 



Fu c’ 


, Fq , Fq , F-^, Fy, F<j> , and F^ become lxnear equatxons 


wxth^constant coeffxcxents xn u c , u s , v c , v s , w c , w s , £ c , £ s , §c> 3s> X, X, 
$ x , and thexr first two txme derxvatxves upon substxtutxon of equa- 

txons (87) . The details of these operatxons are omxtted and only the fxnal 
results are gxven xn the next sectxon. 


7.4 Lxnearxzed Perturbatxon Equatxons of Motxon 

The lxnearxzed perturbatxon equatxons of motion for xnfxnxtesxmal motions 
about equxlxbrxum are expressed xn matrxx form as 

[M 1 + t^HX} + [G 1 + C A + C D ]{X} + [K 1 + K S + K S + + K D ]{X} = 0 (89) 

where 

mass matrxx due to xnertxal forces, symmetrxc 
mass matrxx due to aerodynamxc forces 

gyroscopxc matrix due to inertxal forces, antxsymmetrxc 
damping matrxx due to aerodynamxc forces 


[M 1 ] = 
[G 1 ] = 

[c A ] = 
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[C D ] = damping matrix due to structural damping, diagonal 

[K' 1 '] = stiffness matrix due to inertial forces, symmetric 

[K ] = stiffness matrix due to springs and structural forces, symmetric 

[K®] = stiffness matrix due to gravitational forces, symmetric except as 
specified 

[K^] = stiffness matrix due to structural damping, antisymmetric 

Thus, 1 is the u c inertial term in the u c equation The elements of 
the above matrices are given below. Only nonzero elements are specified, for 
matrices with symmetry or antisymmetry, only the terms on or above the diag- 
onal are given. 

Modifications to the above analysis for different pitch-control geome- 
tries are described m appendices A to C. 
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Recall from section 5 2 that when the aircraft is airborne the parameter 
g is given by Ng 0 where 
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Rows 3 and 4 of are identical to rows 1 and 2, respectively, except 

Pj is replaced by P 2 - 

Rows 5 and 6 of are identical to rows 1 and 2, respectively, except 
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Rows 3 and 4 of C are identical to rows 1 and 2, respectively, except 
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^ 2,10 ” k 1,9 

Rows 3 and 4 of are identical to rows 1 and 2, respectively, except 

that P^ is replaced by P 2 . 

Rows 5 and 6 of are identical to rows 1 and 2, respectively, except 

that P 1 is replaced by P 3 . 

Rows 7 and 8 of are_identical to rows 1 and 2, respectively, except 

that P x is replaced by Q^/jl. 

Rows 9 and 10 of are identical to rows 1 and 2, respectively, except 

that Pj is replaced by QgM. 

Rows 11 and 12 of are identical to rows 1 and 2, respectively, 

except that Pj is replaced by Qq/L 
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k a 
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Rows 15 and 16 of are identical to rows 13 and 14, respectively, 

except that P x and Py are replaced by Q x and Qy, respectively, and the fol- 
lowing terms are added to the elements indicated: 
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The quantities appearing m the elements of the above matrices not defined 
elsewhere in the report are defined as follows: 
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APPENDIX A 


FLEXIBLE TORQUE TUBE 

♦ 

The effect of a pitch-control system similar to that of the Boeing BMR 
(ref. 9) can he easily incorporated into the analysis. This device has the 
effect of placing an additional twisting moment at the tip of the flexbeam. 
This moment is pi'oportional to 8 - where 8 0 is the control input (see 
fig. 4). We assume that the torque tube is sufficiently flexible m bending 
so that no bending moments are applied by the control system at the flexbeam 
tip. With torsion stiffness Kq, the additional moment m the equilibrium 
generalized force expression is 

F 0 = K 0 (0 - 9 0 ) + . . (Al) 

where Kq = Kq/I&q 2 . By changing 0 O , any pitch angle that is desired can be 
obtained from the iterative process described m the text. The stiffness 
matrix is also modified so that 


K 


S 

11,11 


= ^ 12 ,12 " K 66 


+ 


(A2) 


Equations (Al) and (A2) show the modifications necessary to account for a 
flexible torque tube. 
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APPENDIX B 


CANTILEVER PITCH ARM WITH FLEXIBLE PITCH LINK 


A pitch-control system with geometry similar to those of references 7 
and 8 may be easily added to the analysis described in the text. In figure 5, 
a schematic of this configuration is shown, and m figure 12 a more detailed 
schematic of the pitch-link geometry is shown. Here the blade and flexbeam 
are shown with £f = Bf = 6 f = u = v = w= £ = 8 = e = Cb = (3b = 9 b =0 * The 
parameters u l5 v l4 w^ define the position of the swashplate-end of the pitch 
link S with respect to point 0 m the n 1 , n 2 ^, n^ k axis system (there are 
no k’s, indicating that the geometry is the same for all k) . The parameters 
x l 5 Yi* arL d define the pitch- arm/pitch-link junction J with respect to 
point 0 m the n x ^, n y k, n z ^ axis system. Geometric constraints may be for- 
mulated by writing the position vector of P m two different ways One way 
is by the flexbeam tip through the pitch arm: 

P / ° _ \ k k, k. k, k. k /T> , . 

r - (1 + u k )n x + v k n 2 + w fc n 3 + x^ + y x n y + z^ (Bl) 


Another way is by the point S 
P/0 k , k , k . 

r = n 1 n 1 + v 1 n 2 + w-jiij + 


and through the pitch link 

(p 0 + P k ) (sm sin + cos Sin 5 k n 2^ 


+ cos 




(B2) 


Here, all lengths are made dimensionless by Z and the pitch-link length is 
Pq (undeformed) + pk (deformation) The angles Pk and ^k define the orien- 
tation of the pitch link, shown m figure 13. These equations yield geometric 
constraints as follows: 


f x k = ! + u k + XjTJj + yjT^ + ZjT^ - Uj - (p 0 + p k )sm u k sm = C 
f y k = v k + x x T k 2 + Yl T k 2 + z x T k 2 - Vj - (p Q + p k )cos y k sm d k = 0 > (B3) 

f z k = w k + X 1 T 13 + y i T ^3 + 2 1 T 33 " W 1 " (p 0 + P k )cOS \ “ 0 

where f x k is obtained by taking the dot product of (Bl) and (B2) with n x k ; 

fy is obtained by taking the dot product of (Bl) and (B2) with n 2 k , f k is 

obtained by taking the dot product of (Bl) and (B2) with n 3 k 


These geometric constraints form relationships between the geometric vari- 
ables and enable us to eliminate superfluous variables from the analysis. 
Meirovitch (ref. 13) has shown that the generalized forces due to the con- 
straints are 
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k , k 3f x k 3f y k 9 z k _ r „ , ,v 

aq r A x k V k X z k q r V V V ? k’ X 9 k’ P k’ 6 k ,p k 

q r q r q r 

The A x k , A y k , A z k are Lagrange multipliers and the (Xq k ate the generalized 
forces. Thus r 


a Uj__ X x 


a = X 

v k y 


X x z 


% = x i( x y kT n “ x * k ^ 2 ) + yi(^i-^2) + *i(^5i " X x k 4) 

X = - + + “A) + - %Rx + 


a g k = "( X x kc ?k +A y ks C k )( x l T i3 + y l T 23 + z l T 3 3) + X z k [ x i c B k - s 3 k (yi s 0 k +s 
° 9 k = ( yiT31 "" + X y ( y l T 32 “ z 1 ^ 22 ) + X z ^ y 1^33 “ z 1^23^ 

a y k = (P 0 + P k ) sln fi k( A y k Sin y k ” X x k C0S V k) = 0 

a 5 k = (p 0 + p k*[“ cos 6 k( X x k S1T1 u k + ^ C ° S p k) + X z k sin 6 k] = 0 

|{[ — 
a P k = _X x sin Pfc Sln 6 k “ X y cos y k s:m 6 k ” X z C0S 6 k = K p P k 

_ K ** 


lC0 k)] 


Here, = 0 ^ - Kpp k = 0 because the mass and inertia of the link are 

neglected. From ctp = ag = ctp - Kpp k =0, the Lagrange multipliers, which 
are the unknown reactions at the point P, are 

X = "Vk Sin Sin 6k 

1c 

A y = “ K p P k C ° S P k Sin ( B6 > 


k = -K p p k cos 6 k 
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The relations (B6) and (B3) are six relations between the twelve variables 

u k’ v k’ w k’ ? k’ 9 k» y k’ 6 k> Pk’ X x k ’ X y k ’ and x z k ' Thus ’ the blade system 

retains its six degrees of freedom. The scheme for calculating the equilibrium 
deflections needs to be modified slightly. Given uv, v^,, w.j., Xj_, y^-, z^-, and ■■ - - 
Pq. and an estimate of u,” v, w, £, B 5 and 0 we can calculate p, 5, p equi- 
librium components of p^, 6^, p^, respectively, from (B3) and then A x ^, 
and from (B6). Now the total generalized forces including a u , a v , etc., 

from (B5) , are known. Integration along the flexbeam will yield slopes and 
deflections to be included m the minimization scheme, equation (84). The 
iteration may include a variation m w.. , to achieve the desired pitch angle 
or thrust as an extra independent variable m the minimization scheme. 

The perturbation equations are obtained by expressing equations (B3)-(B6) 

m perturbation from and eliminating A x , A x , A y , A y , \ z , A z , u c , u s , 

C S *>C j s c s ^ 

v c , v s , w c , and w s . This operation yields 

[R T MR]{X} + [R T CR]{X} + [R T KR + D]{X} = 0 (B7) 

where 

[M] = -[M 1 + 

[C] = [G 1 + C A + C D ] 

[K] = [K 1 + K S + K g + K*- + K D ] 

{X} — j^P c , p g , 6 C , P c , Pg, 5 C , Cg, 6 C J 8g, 9 C , 0g ,X, Y, , $yj 

0 I 6 0 

0 0 I h 

[I 4 ] = 4x4 identity matrix 

[Ig] = 6x6 identity matrix 
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0 Rn 


0 R, 


0 R" 


0 R" 


o r: 


0 Ro 


l = 1, 2 


o r: 


0 R, 


o rJ 2 


x 1 T 12 + y l T 22 + z 1 T 32 
c C (x 1 T 13 + y l T 23 + z 1 T 33^ 


= 2 1 T 21 - y l T 31 


(Xi T n + YjT 21 + z 1 T 31 ) 

= s ? ( Xl T 13 + yi T 23 + Zl T 33 ) 


z 1 T 22 " y l T 32 


“ - x i c e + s g(yi s e + z i c e> 


z 1 T 23 y l T 33 


= (p„ + p)cos y sm <5 


(pg + p)sxn y cos 6 


= sin y sin 8 


- — (p n + p)sm y sm 8 


= (p n + p)cos y cos 6 


= cos y sin 8 
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R fl - 0 


“ — (p n + p) sin 6 


32 


R 33 


[D] = 


0 

cos 6 


D 1 ol 

,° D J 




D 


11 


0 D 


12 


0 D 


13 


D 


11 


Symmetric 


0 D 


12 


0 D u 


22 


0 D 


23 


22 


D 


0 D 

1 

33 


>i. ■ 


D 


12 


(P 0 + P) ( K pP + cos <$) 
0 


D 13 " 0 


D 22 = ( P 0 + P)K P P 


D 2 3 0 


D 


1 _ 

33 


= S 


D U = " A x R li + Vll 


D ?2 


D 2 

u 13 


^x 1 *^ + ^y R l2 

~^x R 23 + ^y R 13 


23 

0 

j 1 

J 33 


(B8) 

con- 

tin- 

ued 




) 2 = 
J 22 


= -(A x c c + A y s ? )Rl 2 + A z (c ? Rl 2 + s r Rl,) 


3 C 22 ; 
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o*jto cj'll'O 


D 2 3 ^x c ? + V S ?^ R 3 3 + 3 "[" S £?2’3^ 

D 33 “ ^X^1 T 21 + Z 1 T 31^ + ^y^ y l T 22 + Z 1 T 32^ + *2^1^23 t z lT33^ 



Ay + A. 

A 2 + A 
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APPENDIX C 


SNUBBER AND PITCH ARM WITH FLEXIBLE PITCH LINK AND SNUBBER CONNECTION 


For the case of a pitch arm and snubber as in the Sikorsky UTTAS tail 
rotor (ref. 6), the geometry may be treated as two flexible pitch links with 
two pitch arms and two different sets of geometric variables* u^ , , w^ , x^ , 

y j , Zj, and Pq^J and u 2 , v 2 , w 2 , x 2 , y 2 , z 2 , and Pq 2 ‘ The two dimensionless 
flexibilities are and K 2 and the deflections, analogous to appendix "B, 

are and p 2 ^, d 2 ^-, p 2 ^. A typical configuration of this type is 

shown m figure 6. 


The mathematical development of the generalized forces and the iterative 
equilibrium solution is identical with that of appendix B except that there 
is a respective contribution from both pitch arms. 

,-k i . . _k , _k , _k , , k. k . k _ 

f Xl = 1 + u k _ U 1 +x l T ll +y l T 21 +z 1 T 31 “ (Poj+Pl )sin P 1 Sln 6 i = 0 


: x 1 + u k U 2 +X 2 T 11 + y 2 T 21 + z 2 T 31 


etc 


( , k. k „ k n 

(Pq 2 + P2 ) sln w 2 Sln 5 2 = 0 


_,k , .k 

a u, A xi + A X2 


etc 


a 


« k ■ + yi»y 1 T 21 - 4,4) +z d4 1 T 31 - 

+ V4 2 T U - 4 2 4 2 > +y 2 «y 2 T 21 - 4 2 T 22> + z 2<4 2 T 31 - 4 2 T 3 2 ) 


etc 


•, k _ k »k-k 

A y = -sm y_ sin S, K p 

i i i l i 


l = 1, 2 


etc. 


(Cl) 


Therefore, the equilibrium scheme is altered m a way similar to that of 
appendix B except that both pitch arms contribute. One pitch link will, of 
course, change the pitch angle where the other one remains fixed as a snubber. 


The perturbation equations are given by 

[M] {x} + [C]{x> + [K + Q]{x} = 0 


(C2) 


74 


ORIGINAL PAGE1S 
OF FQGR QUALff*i 



where 



where the subscript i on the matrix [D^] refers to subscripting all the 
parameters in [D 2 ] sub i (x x , y 15 Zj_, A. Xi , etc.) and similarly for 

the following 


[SJ » -[R 1 ] 1 i = 1, 2 

[EJ = [A x T C ^ 1 Aj i = 1 , 2 

[Aj - [R 2 ] 1 i = 1, 2 

[C 1 ] = [D 1 ] 1 i,= 1, 2 


(C4) 


See appendix B for the definitions of the matrices on the nght-han4 side of 
equation (C4) . 
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APPENDIX D 


STRUCTURAL DAMPING COEFFICIENTS 


In section 5.5, the blade and body structural damping coefficients are 
used m the analysis. We first consider the rotor blade for Q = 0 m vacuo 
It is desirable to specify a structural damping ratio isolated blade chordwise 
motion. The single blade motion is given by 


m£ 2 

mJlx^ 


[♦1 


0 

0 


fvl 


12B 

-6B 


(i 



* 


• + 




• 

' + 



■ 


inJlx^ 

1 


[~ ? J 


0 

C c_ 




-6B 

4B 




(Dl) 


where c^ is unknown. We wish to choose c^ so that the motion is critically 
damped and then multiply by a dimensionless number given by the struc- 

tural damping ratio. For critically damped motion we require that the eigen- 
values of (Dl) be a pair of complex conjugates -a ± ico and a pair of negative 
real numbers -a, -a. The characteristic equation for such roots is 

s 1 * + 2 (a + n) s 3 + [ (a + a) 2 + 2aa + to 2 ] s 2 + 2 [an (a + cr) + ato 2 ] s + a 2 (a 2 + to 2 ) = 0 (D2) 


The characteristic equation for equation (Dl) is 

+ p 3 Cj.s 3 4- p 2 s 2 + PjC^s + p Q = 0 


where 


P 0 = PjB 

12B 

Pi = — 

m£ 2 (l - mx^ 2 ) 

/. -- m& 2 \ 

p 2 = p i v 1 + m£x b + ~ir) 


p 3 



) 


(D3) 


(D4) 


When the coefficients of equations (D2) and (D3) are equated, it can be shown 
that 


3a 4 + p 2 a 2 - p Q 
2a 3 p 3 


(D5) 
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where a 2 xs the sxngle real, posxtxve root of 

^ 3p l 


a 6 + a 4 * 


P3 ' + " 


i / P l P 2 

\ P 3 


3Pr 


P 1 P 0 


= 0 


(D6) 


It xs dxfficult to achxeve an accurate numerxcal solutxon, for small £, to 
equatxons (D5) and (D6). An asymptotxc expansxon, however, presents no prob- 
lem numerxcally and yxelds a sxmple formula accurate to second order xn l’ 


= v'B 2 + mlx b + ^1 - ix b 2 ) + 0(£ 3 ) 


(B7) 


which gxves accuracy wxthin 1% for £ < 0.27. 
The structural dampxng value xs thus 


c 5 = 


-~2 / 4mx, 

2 + S£x b + (l 3 


where typxcally 0.005 < p^ < 0 03 

For the body motxon, we choose structural dampxng coeffxcxents as xf 
X, Y, $ x , and Oy motxons were uncoupled* 


” 2p X 3 , 1 3 M 1 3 , 1 : 


" 2p y 1/k 3 4) i 4 m[ 4> i4 
2t1 $ v V K 15,15 M 15,15 


c $ = 2n$ 
y *y 


1^ K 16 > 16^16 > 


16 


Here Py, q<| , and p^ are the body structural dampxng coeffxcxents, 

x y 

typxcally 0.005 < p x , p y , p $ , p $ <0.05 
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Figure 1 - Rotorcraft fuselage model. 
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Figure 


2 - Landing gear geometry. 
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Figure 3.- Rotor blade configuration. 



Figure 4.- Flexible torque tube. Case II. 
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Figure 5.- Cantilever pitch arm, Case III. 


FLEXBEAM 



Figure 6.- Cantilever pitch arm with snubber. Case IV. 
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BLADE-FIXED 

AXES 


Figure 7.- Rotating rotor blade coordinate systems. 





Figure 8.- Sequence of angles describing blade orientation. 
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Figure 11.- Longitudinal plane of body-spring configuration. 
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Figure 12.- Detailed schematic pitch link configuration. 
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